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Due to the orthorhombic distortion of the lattice, the electronic hopping integrals along the a 
and b diagonals, the orthorhombic directions, are slightly different. We calculate their difference in 
the LDA and find t'^ ~ t'f, ~ 8 meV. We argue that electron correlations in the insulating phase of 
La2-a:Sra:Cu04, i. c. at doping x < 0.055, dramatically enhance the {t'a — tf,)-splitting between the a- 
and 6-hole valleys. In particular, we predict that the intensity of both angle-resolved photoemission 
and of optical absorption is very different for the a and b nodal points. 



INTRODUCTION 



The magnetic state of La2-2:Sra;Cu04 (LSCO) changes 
tremendously with Sr doping. The three-dimensional an- 
tiferromagnetic Neel order identified [l[ below 325 K in 
the parent compound disappears at doping x ~ 0.02 
and gives way to the so-called spin-glass phase which 
extends up to a; w 0.055. In both, the Neel and the 
spin-glass phase, the system essentially behaves as an 
Anderson insulator and exhibits only hopping conduc- 
tivity [l|, Q. Superconductivity then sets in for dop- 
iiig ^ ^ 0.055, see Rcf. [H One of the most intriguing 
properties of LSCO is the static incommensurate mag- 
netic order observed at low temperature in elastic neu- 
tron scattering experiments. This order manifests itself 
as a scattering peak shifted with respect to the antifer- 
romagnetic position. Very importantly, the incommen- 
surate order is a generic feature of LSCO. According to 
experiments in the Neel phase, the incommensurability 
is almost doping independent and directed along the or- 
thorhombic b axis In the spin-glass phase, the shift 
is also directed along the b axis, but scales linearly with 
doping 0, 0, Sj- Finally, in the undcrdopcd supercon- 
ducting region (0.055 ^ x < 0.12), the shift still scales 
linearly with doping, but it is directed along the crystal 
axes of the tetragonal lattice In the present work 
we discuss only the insulating phase, x < 0.055. It is 
clear that pinning of the diagonal spin structure is due 
to the orthorhombic distortion of the crystal. A mech- 
anism for pinning of the diagonal spin structure to the 
orthorhombic b axis was suggested in Refs. The 
mechanism has four components: 

1) Due to strong antifcrromagnetic correlations, the 
minima of dispersion of a mobile hole are at points 
(±7r/2, ±7r/2) of the Brillouin zone, so the system can, 
to some extent, be considered as a two valley semicon- 
ductor. 

2) At low temperature, each hole is trapped in a 
hydrogen-like bound state near the corresponding Sr ion, 
the binding energy is about 10 meV and the radius of the 
bound state is about 10 A. 



3) Due to the orthorhombic distortion, the diagonal hop- 
ping matrix elements t'^ and are slightly different, and 
this makes the b valley, (— 7r/2, 7r/2), deeper than the a 
valley, (7r/2,7r/2). So all the hydrogen-like bound states 
are built with holes from the 6- valley. 

4) Each hydrogen like bound state creates a spiral distor- 
tion of the spin background and the distortion is observed 
in neutron scattering. So the state at 0.02 < x < 0.055 
is not a spin glass, it is a disordered spin spiral. 

In the present paper we calculate accurately the diago- 
nal hopping matrix elements and show that the difference 
is tjj — « 8 meV. We also discuss the implications of 
the described physics for the p^cP p^(P optical ab- 
sorption and for angle-resolved photoemission (ARPES). 
The optical absorption probes the charge distribution, 
and we predict that the absorption vanishes for polariza- 
tion along the orthorhombic a-direction. ARPES probes 
the spin spiral, and we predict very different spectra in 
a and b nodal points. 



ANISOTROPY OF t' DUE TO ORTHORHOMBIC 
DISTORTION AND TILTING 

The 2D < — J model was suggested two decades ago 
to describe the essential low-energy physics of high-Tc 
11, 12, 13|. In its extended version, this model 



cuprates 

includes additional hopping matrix elements t' and t" 
to respectively 2nd and 3rd-nearest Cu neighbors. The 
Hamiltonian of the t — t' — t" — J model on the square 
Cu lattice has the form: 

{ij)a {ij")^ 



(1) 



Here, is the creation operator for an electron with spin 
a (a at site i of the square lattice, (ij) indicates 

1st-, {ij') 2nd-, and {ij") 3rd-nearest neighbor sites. The 
spin operator is Si = ^cl^cTa^Ci/j, and Ui = J^a^a^^io- 
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with {rii) = 1 — X being the number density operator. In 
addition to the Hamiltonian ^ there is the constraint 
of no double occupancy, which accounts for strong elec- 
tron correlations. The values of the parameters of the 
Hamiltonian ([T]) for LSCO are known from neutron scat- 
tering [1|, Raman spectroscopy 14 1 , and ab-initio calcu- 
lations [l5| to be: 



J w 140 meV, 
t' « 70meV, 



450 mcV, 
35 meV. 



(2) 



Note that the signs of the hopping terms in the electron 
Hamiltonian (JT]) have been chosen in such a way that, 
for a (1^2 _y2 orbital, the hopping matrix elements are 
positive. 

The dispersion of the hole dressed by magnetic quan- 
tum fluctuations has minima at the nodal points qo = 
(±7r/2, ±7r/2) and is practically isotropic in the vicinity 
of each d: 



e(q) w e(qo) + ^/3(q-qo)' 
P « 2J« 260meV . 



(3) 



We set the lattice spacing to unity, 3.81 A —> 1. The effec- 
tive mass corresponding to the quadratic dispersion ([3]) 
is approximately twice the electron mass. 




FIG. 1: Orthorhombic deformation of the square lattice 

In the low-temperature orthorhombic (LTO) phase, 
the square Cu lattice is slightly deformed as indicated in 
Fig- El The angles, 7r/2±(y5, between the edges are obtuse 
or acute, such that the (11) and (1 — 1) diagonals, the or- 
thorhombic translations, have slightly different lengths, 
respectively a = a/2 (1 - (/3/2) and 6 = \/2 (1 -t- <^/2) . On 
top of this comes an alternating tilt of the oxygen oc- 
tahedra around the a axis by an angle ±9, whereby the 
oxygens in the layer buckle out of the Cu plane forming 
a [— 7r/2, 7r/2]-wave in the 6-direction [H, [l^. The values 
of the orthorhombic and tilting deformations. 



(^ = 0.009 (0.5°) 



and 



e = 0.05 (3°) , (4) 



are so tiny that they hardly influence J, t and t" , but 
they do make the diagonal hopping significantly different 
in the a and b directions, as we show below. Including 



t'^ t'^ in the first, one-electron part of the Hamiltonian 
([1]) yields the following 2D bandstructure: 

e (k) = — 2t (cos kx + cos ky) — 2t" (cos 2kx + cos 2ky) 
+2t'^ cos {hj, + ky) + 2t'i, cos {k^ - ky) . (5) 

Note that e denotes electron energies and e hole energies. 
The splitting of the bare valleys is: 



e(-7r/2,7r/2)-£(7r/2,7r/2) 



£b 



(6) 



4 (t; - t',) EE 46t'. 



For comparison, we show in the left-hand side of Fig. [2] 
the bandstructure calculated from density-functional the- 
ory (LDA) along the lines connecting the high-symmetry 
points r (0,0,0), R (7r,0,7r/2), S (7r/2, 7r/2, 7r/2) , Y 
(0, 0, 7r/2) , and B (-7r/2, tt/2, tt/2) of the 3D oHhorhom- 
bic (Bmab) Brillouin zone (BZ). For simplicity, we have 
given the coordinates in terms of the nearly tetragonal 
reciprocal Cu-lattice translations, kz = tt/2 is halfway 
to the Brillouin-zone boundary where the influence of 
inter-layer hopping is minimal. We see that no gapping 
is caused by the [— 7r/2,7r/2] tilting wave. As a conse- 
quence, the LDA bandstructure may be folded out to 
the nearly tetragonal BZ where expression ([5|) adequately 
represents the dispersion of the LDA conduction band, 
indicated by heavy lining, near half-filling. On the righ- 
hand side of Fig.[2]we show a blow-up of the orthorhombic 
LDA bands near the nodal points, S (7r/2, ■7r/2, 7r/2) and 
B (7r/2, — 7r/2, ■jt/2) , along the respective nodal direction, 
a and b. Finally we see that two different computational 
techniques, LMTO and LAPW, give essentially the same 
result, namely 



St' w 8meV, i.e. dt'/t' w 11%. 



(7) 



This number is surprisingly large: Had the Cu o?^2_y2 
conduction-band Wannier-like orbital been a simple, 
canonical orbital |2l| , the hopping between two such or- 
bitals would have decreased as their distance to the power 
— {I + I' + 1), which for two d orbitals is —5, and, hence, 
St' /t' ^ 5lp ~ 4.5%. This, however, neglects that the 
hopping is almost exclusively via the O^: Px and OyPy 
orbitals, as is amply demonstrated by an LDA calcula- 
tion in which we took the tilting wave to run along a 
instead of along b. The result was a four times reduc- 
tion of St' compared to ([7])! We therefore conclude that 
tilting is almost as important as orthorhombicity for the 
hopping anisotropy. 

To understand why, we turn to the simplest model 
showing an effect of tilting: Emery's 3-band model [l^ 
which keeps merely the tpd hopping between Cu d^^-y^ 
and its nearest 0^ Px and Oj, Py orbitals, as well as the tpp 
hopping between the nearest Oa; Px and OyPy neighbors. 
It is easy to see that in this model. 



t' 



. tpdtpptp^i 

(e_F - £p) 



so that 



St' St 



^ ^ ^ (8) 
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FIG. 2: LDA bandstructure of the LTO phase shown in 
the orthorhombic BZ (see text), (a) NMTO downfolded Cu 
d^2_y2 - like conduction bands shown in thick Unes [3]. (b) 
Blow-up of conduction bands near the nodal points S and B, 
and along the respective nodal lines, a and 6. (c) same as 
(b), but calculated with the full-potential Linear Augmented 
Plane Wave Method (LAPW) [MJnstead of the Linear Muf- 
fin Tin Orbital Method (LMTO) ^ 



because tpd and the distance of the O p level, Sp, below 
the half-fiUing level, ep, are not influenced by the defor- 
mations. In the deformed structure, each plane oxygen 
remains midway between its two nearest Cu neighbors, 
but is slightly above or below the Cu plane. Moreover, 
the Op orbital is parallel to the Cu-Cu line, so that the 
two Op orbitals via which the t'^^^ hopping takes place 
are at the acute/obtuse angle tt/2 ^ (p. The main effect 
of the deformation is, however, not this misalignment, 
but simply that the distance between the oxygens at the 
acute angle is a/2 and that between the oxygens at the 
oblate angle is b/2/ cosd. With I = I' ~ 1, this gives: 
dt'/t' = 3 (v? + 6*2) = 3 (0.90 + 0.25) % « 3.5%, which, 
due to the slower decay of canonical p orbitals, is even 
smaller than the previous dd estimate. Also the tilting 
contribution is too small. Inclusion of the orbital mis- 
alignments in the canonical approximation [2H merely 
changes the result to: 6t' /t' = ^ip+bO'^ w (2.1 + 1.3) % = 
3.4%. 

Another simple model is the axial-orbital model [l^ 
which adds to the three orbitals of the Emery model a 
Cu-centered, axial orbital, but keeps only hops between 
nearest neighbors: from Op to Cu dx^-y^ (tpd) and to 
the axial orbital {tsp) ■ The axial orbital is a hybrid be- 
tween Cu4s, Cu 3(^322 _]^, apical Opz and axial cation 
orbitals, and its energy, (> ep) is the material depen- 
dent electronic parameter in the LDA. Since tpp proceeds 
via the axial orbital in this model, and therefore equals 
^Isp/ i^s — £f) , it cannot depend on ipl For the same rea- 
son, the hopping from Ox Px to O^, Px over the distance 1 
has the same value, tpp, and this is what leads to t" ~ t'/2 
when downfolding to the 1-band model [isj . 



However, recent first-principles 3-band Hamiltonians 
formed by numerical downfolding of the LDA Hilbcrt 
space [18| do not support this; they yield a Px-Px hop, 
which can only be understood by the presence of a large 
material-independent contribution, — t|pp/ (e4p — £f) , 
from hopping via Cu Apx [l^]. To the Px-Py hop, there 
can be no such contribution, — unless there is orthorhom- 
bic distortion. Also buckling-induced anisotropy may be 
caused by hopping via Cu ipz- 

To see whether coupling via excited Cu 4p degrees 
of freedom can be the reason for the surprisingly large 
anisotropy ([7]) found by the LDA calculations, we now 
add these degrees of freedom to the axial model and, as 
usual, include hops only between nearest neighbors. For 
the directional dependences we use the canonical approx- 
imation and find: 



5t^ 



St 



pp 
tpp 



= (2^ + 4.56*2 



4pp 



(1.8 + 1.1)% X 



sp 

2.6 
2^ 



£4p 
2 



35 
IT 



12%. 



The superb agreement with the LDA result ([7]) may be 
fortuitous, but also the relative contribution from or- 
thorhombicity and tilt agrees. So we believe that the 
surprisingly large hopping anisotropy to be explained by 
this simple expression. It may be noted that going to 
materials with higher Tcmax, decreases and with it the 
hopping anisotropy. Finally it should be noted that for 
simplicity we have taken £4^ = £4^ (= £42:/y) , although 
£4z is an axial orbital and, hence, material dependent. 

The t' describes hopping within the same magnetic 
sublattice. Therefore, to account for many electron corre- 
lations in the anisotropy of the single hole dispersion one 
only needs to know the quasiparticle residue Z (q) . This 
leads to anisotropic correction to the hole dispersion Q 



Se (q) = 2Z (q) 5t' sinqx sinqy 



(9) 



The correction to the dispersion vanishes at antinodal 
points (0,7r), (tt, 0) and it is maximum at nodal points, 
(±7r/2, ±7r/2), where Z w 0.3. The energy difference 
between the nodal points is 



£6 



AZSt' - lOmeV 



(10) 



Thus, at low temperature all hydrogen-like bound states 
are formed from the 6- valley holes, the a-valley is empty. 
We would like to stress that this is true only in the insu- 
lating phase. In the superconducting phase which arises 
after percolation of the bound states, x > 0.055, both 
valleys are populated. 

The energy difference ^TU\\ is the single hole effect. It 
accounts for spin quantum fluctuations, but it assumes 
usual Neel order. There is another collective contribu- 
tion to the energy difference — £{,. The collective con- 
tribution is due to the spiral spin ordering established at 
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0.02 < X < 0.055. The collective contribution is consid- 
ered below in the ARPES section. Certainly, in the end 
the collective contribution is also driven by the asymme- 
try PT?|) because the asymmetry's causes depopulation of 
the a-valley. 



IMPLICATIONS FOR OPTICAL TRANSITIONS 

We consider the optical transitions involving removing 
a d electron from a central placket of Fig. [3] resulting in 
a d^-state and transferring it to a linear combination of 
Zhang Rice (ZR) singlet states on neighboring plackets. 
Note that the transition staying within the same plaque 
is not allowed optically. This particular transition is of 
course only allowed if there is a hole in the ZR state on 
neighboring plackets and therefore is only present in a 
hole doped material. In the undopcd material the cor- 
responding transition would have to be to neighboring 
d states resulting in a linear combination of d^'^ config- 
urations there. This transition would be at an energy 
of about the charge transfer gap higher than the tran- 
sitions to ZR singlet states and is easily distinguishable 
from those to the ZR states. These transitions involve 
the full d^ multiplet structure and their energies can be 
obtained from the calculations described in Ref. [i^ as 
well as from early resonant photoemission experiments 



24| 



which locate the d^ multiplets in these materials. 
Looking at Fig. 1 of Ref. [23[ we see that the d^ states 
occur in the energy range between 8 and 15 eV below the 
ZR singlet state. The large energy spread is due to the 
multiplet structure in the d^ configurations. 

We now describe why the polarization dependence of 
these transitions will be very sensitive to where in mo- 
mentum space the ZR states are situated and therefore 
the polarization dependence represents a test of the hy- 
pothesis described above. For this discussion the or- 
thorhombic distortion is no longer important except for 
the proposal that it causes the doped hole states to be 
concentrated at the b minimum with a wave function 
given by 



V'b(r) = X(r)exp |i|a; - i|y| , 



(11) 



where x(r) oc e~^'' is a relatively smooth wave function 
of the bound state, k « 0.4. The smooth wave function 
is not important for optical absorption. Only the fast 
phase factor exp {i^x — is important. For a bound 
state based on a hole from the a-minimum the function 
x(r) is the same while the phase factor exp [i^x + i^y} 
is different. 

To describe the above transition we clearly need 
to go beyond a. t — J or Hubbard model taking 
into account explicitly the charge transfer nature of 
the gap as in the Zaancn-Sawatzky-AUcn classification 



scheme [25| . In Fig. [3] we display the Cu02 plane struc- 
ture and the orbitals considered. The hole phase factor 
exp {i^x — i^y} = ±1, ±« is shown in red near the cor- 
responding Cu ion. The transition to a zero momentum 
ZR state is parity forbidden. However the hole resides in 
a ZR state with nonzero momentum where transition is 
optically allowed due to interference of ZR states centered 
on different Cu sites. Let us consider the case when the 




+ 1 



+1 



-1 



FIG. 3; (Color online) Cu02 plane. Open circles denote Cu 
ions and blue circles denote Oxygen ions. Cu ions in the 
cluster are enumerated by magenta numbers 1-9 and Oxygen 
ions are enumerated by four magenta letters A,B,C,D. The 
phase factor, ±1, ±i, is shown in red near the corresponding 
Cu ion. 

final hole resides at the first Cu site, \ final) ~ \di). Here 
for short notations we denote d^ state as d. The transi- 
tion p — !■ c?i can occur when the initial p-holc resides on 
Oxygens A,B,C,D and comes from the ZR states centered 
at Cu sites 2,4,6,8. Therefore the transition amplitude is 

A = (di \ExX + Eyy \ - ipA + ipc - Wb + ipo) , (12) 

where E is the electric field of the photon. Let us denote 
by D the dipole matrix element 

{di\x\pc) = X> = J ipdiXippcdV . (13) 

The symmetry relations between nonzero matrix ele- 
ments follow from Fig. [3] 

{di\x\pA) = -T) , 
{di\y\pB)=V , 

{d,\y\pD) ^ -V . (14) 

Hence, we find from 

A = 2i{E, - Ey)V . (15) 

Thus we conclude that there is an absorption if the wave 
is polarized along the orthorhombic 5-axis, E oc (1,-1) 
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and there is no absorption if the wave is polarized along 
the orthorhombic a-axis, E oc (1,1)- Basically this is 
the only possible correlation one can write kinematically, 
I (X {E ■ fc)^, where k is momentum of the hole. So, 
the answer is obvious even without a calculation. The 
optical absorption we have discussed is proportional to 
doping X. We stress that the prediction for the low tem- 
perature absorption asymmetry is ab ^ Ua and this is 
equally applicable to the Neel and "spin-glass" phases of 
LSCO. The optical asymmetry is different from that in 
the dc conductivity where both experimentally and 
theoretically the asymmetry at low temperature is 
not that large, ^50%, and of the opposite sign, ai, < Ua. 

IMPLICATIONS FOR ARPES 

We consider photoemission from a and h nodal points. 
There is a difference in energy that is given by Eq. ()10p . 
This is already an interesting effect. However, there is a 
much bigger effect that is due to the spin spiral in the 
"spin-glass" phase. The pitch of the spiral directed along 
the 6-axis is [1, 

Q = ^(l,-1), (16) 

Ps 

where ps ~ 0.18J is spin stiffness, g Zt 0.7J is 
the hole-spin- wave coupling constant, and x is doping. 
Eq. (jl6p agrees very well with neutron scattering data. 
A hole with momentum q interacts with the spiral. The 
interaction splits the hole dispersion in two branches with 
the following energy shift [§, [l^ 

Ae (q) = iplQj; sing^ + Qj, singj^l . (17) 

Thus the hole dispersion near the "a" nodal point, q w 
(7r/2, 7r/2), is practically not influenced by the spiral, 
Aea = 0. On the other hand the hole dispersion near 
the "6" nodal point, q w (7r/2, — 7r/2), is changed as 

Ae, = ±^ . (18) 

Ps 

Accounting this correction together with pI7|) we find 
eb - ea = -4Z(5i' ± 10±770x meV . (19) 

Ps 

Note that the second (collective) contribution of this dif- 
ference scales linearly with doping and it is pretty large, 
it is --30meV at x = 0.04. 

Thus the present picture predicts that the lowest 
branch of the dispersion with energy approximately equal 
to chemical potential is at the nodal 6-point. The dis- 
persion at the nodal a-point is by ~40meV higher (we 
present the estimate for x = 0.04). Finally, there is an- 
other branch of the dispersion at the nodal 6-point that 



is by ^'GOmeV above the chemical potential, this branch 
can be pretty broad. The present consideration is ap- 
plicable to the doping interval 0.02 < a; < 0.055 where 
the disordered spiral is established. In the Neel phase, 
X < 0.02, the collective contribution in is suppressed. 

CONCLUSIONS 

We have demonstrated that the 0.9% orthorhombic 
distortion of the Cu lattice causes 7% difference in diag- 
onal hopping matrix elements, and that the 3° tilting of 
the oxygen octahedra causes an additional 4% difference. 
This 11% difference together with the effect of strong 
magnetic fluctuations (small hole pockets) and together 
with localization of holes due to Coulomb trapping by 
Sr ions leads to depopulation of the (7r/2,7r/2) pocket in 
the insulating phase, x < 0.055. As a result the optical 
transition p^dP p^(P is allowed only if the electric field 
is polarized along the orthorhombic 6-axis. 

Another prediction that is related to the spin spiral 
structure is the asymmetry of ARPES spectra: the hole 
dispersion at the 6-nodal point is close to the chemical 
potential, while at the a-nodal point it is by 25-50 meV 
higher depending on doping. 
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